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1 . Abstract 

A model reduction set hod for discrete bilinear systems la developed which Batches q sets of 
Volterra and covariance paranetera. These parsnetera are shown to represent both deterministic 
and stochastic attributes of the discrete bilinear sytten. A reduced order model which matches 
these q sets of parameters la defined to be a q-Volterra cov ariance equivalent .realization (q- 
Voltcrra COVER). An algorithm la presented which construct* a class of q-Volterra COVERS 
parameterized by solutions to a Hermltlan, quadratic, suitrlx equation. The algorithm la applied 
to a bilinear model of a robot manipulator. 




2. Inroductlon 


While model reduction of linear systems has been extensively researched over the past few years, little work 
haa been done In the area of model reduction for nonlinear systems. One class of nonlinear systems which Is 
especially appealing are bilinear systems ((ll-(A]>. Bilinear systems are linear In the state variables, linear 
In the control variables, but nonlinear In the state and control. One reason that this claas Is of Interest Is 
that nonlinear systems which are linear In the control variables can be accurately approximated by bilinear 
models ([5], [6]). Bilinear approximations will In general have a higher order than the original nonlinear system 
and effective means for reducing bilinear models are needed. 

Host approaches to model reduction of linear systems have strived to preserve or approximate a certain 
characterizing property of the full order model. For deterministic systems this property Is typically the 
Impulse response sequence or the system Hankel matrix (e.g., [7]— [101). Hod el reduction of linear stochastic 
systems usually Involves the output covariance sequence or the corresponding Hankel matrix (e.g., Ill] and [12]). 
A model reduction technique which considers both deterministic and stochastic properties has also been developed 
<[131— (151) and the resulting reduced order models have been called q Harkov COVERS ( cova riance equivalent 
realizations) . The model reduction problem for discrete bilinear systems has recently received some attention. 
Hsu et al . (16) develop a method for deterministic, discrete, bilinear systems using a generalized Hankel matrix. 
Deaal has proposed an approach to stochastic model reduction based on his realization theory ([17]). 

In this paper we develop a model reduction algorithm analogous to the q Harkov covariance equivalent 
realization approach for linear systems. The algorithm produces a class of reduced order models which exactly 
match a specified number of deterministic and stochastic parameters. This class of reduced order models Is 
parameterized by the solutions to a Hermltlan, quadratic, matrix equation. Section 3 presents the deterministic 
and stochastic attributes of a bilinear system which we will preserve In our method and defines a q- Volterra 
covarianc e equivalent realization . Next, In section 4 the model reduction algorithm Is outlined. In section 5 a 
parameterization of reduced order models which match q Volterra parameters sod q covariance parameters Is 
formulated. Section 6 contains an application of Che proposed algorithm to a two degree of freedom robot 
manipulator. The final sections are our concluding remarks, acknowledgements and references. 


3 . q-Volterra Covariance Equivalent Realization s 

Consider the time Invariant discrete bilinear system 

n 

u 

x(k+l) - Ax(k) + Z (N x(k) + b.)u.(k) 

1-1 1 11 

y(k) - Cx(k) (1) 

where A and N, , 1-1,..., n are n *n matrices, b, , 1-1,..., n are n »1 matrices and C Is an n »n matrix. The 
1 ’ u xx ’ 1 u x yx 

state vector x(.) Is “ I , the Inputs Uj(.), 1-1, ...,n u are scalar, zero mean. Independent Guassian white noise 
processes with Eu^(J)u^(k) - 4 ^ and for j » k, Ex(k)u^(j) • 0. The output y(.) is an n^«l, zero mean, 

stationary stochastic process. We assume that the bilinear system driven by unit Intensity Guassian white noise 
Is stable In the sense that the state covariance 
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X A llmEx(k)x*(k) > 0 (2) 

k*“ 

Is finite. It esn bt shown that for these Input processes the state covariance will satisfy the bilinear Liapunov 
equation 

n 

fk O A A i 

X - AXA + I H XH. ♦ BB , B « l b, ... \ J . (3) 

1-1 1 1 1 J u 

U. also assume chat Chara ara no redundant lnputi or outputa (B has linearly Independent coluana and C has 
linearly Independent rows). 


Xubertl at al. (3) define the product • for an s»l vector a and an r*l vector b and the product (] for - an 
n«nr matrix L and an n*a matrix M by 

. L 0 H - { Lj ... L r J □ M - ( LjM ... L r M J . 

They also eeteblleh the following Identity, 

L( (Me » b) ] - (L □ M)(a * b) . (*} 


a * b 


ab. 


•b_ 


With these definitions (1) becomes 

x(k+l) - Ax(k) + N[x(k) * u(k)l + Bu(k) 

y(k) - Cx(k) (5) 


and (3) may be expressed as 


X - AXA* ♦ (MQX)H* + BB* , N $ I N, ... N 1 . 

1 a 


(6) 


The xero initial state response of the bilinear system (5) Is sn Infinite Voltarra series [A] . This series 
In regular form Is found to be 


y(k) 


^ 1 "1 ^ k-lj i 1-1 

£ CA 1 Bu(k-lj ) + r I CA 2 NQA 1 Blu(k-l l -l 2 )*u(k-l 2 )l + 


V 1 


V‘ l !-‘ 


k k-l 3 k-l 2 -l 3 t t .j l-l 

Z Z Z CA nCJA N0a * B[u(k-l 1 -i 2 -l 3 )e u (k-l 2 -l 3 )*u(k-l 3 )l ♦ ... 

1 3 -1 1 2 -1 1,-1 


where Identity (A) has been used repeatedly. The matrix valued function In each of the summations It called a 
Volterra kernel, the J 1 * 1 Volterra kernel In regular form Is then 


1-1 1 .-1 i.-l 

hjClj.lj-j l u i CA J NClA J 1 N...NQA 1 B 


(7) 


where 1 >1, n-l,...,j and the matrix N occurs exactly 1-1 times. The step response, u(k)-l for all kX), la 

b n 


y(k) 


k-i. 


J-i ij-i ij.,-1 


• • «“ij 


V 1 


WVl-i’ 


•*i )l j 


where 1 la a column vector of ones with a elements. We shall call the coefficients in the step response the 

Volterra parameters of the bilinear system (5). The Volterra parameters of the Volterra kernel are n «n^ 
c . y u 

matrices. We define the set of q order Volterra parameters as those coefficients In which the aatrlces A and N 
occur a total of q times. For example, 

V 2 - ( CA 2 B , CANDB , cnDab , cnQnQb ) . 

We now see that Che step response Is completely characterized by Che sets of Volterra parameters. We also observe 
that for each k a new set of Volterra parameters effects this response. That is. If a reduced order model matches 
Che first q sets of Volterra parameters of the full order model then It will also match the step response for 
k-0,1 q+1. 


In addition to Volterra parameters we are concerned with a covariance sequence for the bilinear system. 
Desai (17) an! Frazho (18) utilize a covariance sequence which Includes both second moments of the output and 
higher moments between the output and Input processes In their bilinear stochastic realization theories. We also 
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uaa this typa of aaquence, in particular tha acquanca of concern Is 


K,j<0) - Ey(k)y*(k) - CXC* 

R,<1) - Ey(k+l)y*(k) - CAXC* 

Rj(0,0) - Ey(k*-l)[y(k) • u(k)]* - CNQXC* 

R 2 (2) - Ey(k+2)y*(k) - CA 2 XC* 

R 2 (0,1) - Ey(ks-2)|y(k) • u(k,l))* - CNDAXC* 
R 2 (1,0) • Ey(k+2)[y(k) • u(k)]* - CANQXC* 
R 2 (0,0,0) - Ey(k+2)[y(k) • u(k) • u(k+l)]* - CNDdOXC* 


where the subacrlpt Indicates the total number of occurrences of A and N, and the Integer* In parenthesis 
represent the powers of A from left to right. A typical element of the sequence Is then 

R j-l+l J +l J _ l 4-...+l 1 <1 j ,1 J-l 1 1 > 

^ Ey<k+j-l+l J +l J1 + ...+l 1 )|y(k)*u(k+l 1 )*u(k+l-H 2 +l l )*...*u(k+J-2+l j _ i + ...+i l )]* 

i. i... i, . 

- CA j NQA j N...NOA XC . (8) 

As with the Volterra parameters we shall define the set of q c>1 order covariance parameters, K , as thoss 

q 

covariances In which the matrices A and N occur a total of q times, that Is the set of second order covariance 
parameters Is 

r 2 - { ca 2 xc* , cnOaxc* , canQxc* , CNDNDXC* } . 

These sets of covariance parameters completely characterise the stochastic bilinear system. It Is worth noting 
that If a reduced order model matches tha first q sets of covariance parameters of the full order model It will 
also match exactly the mean square value of the output and all output and Input correlations up to q steps Is 
time. 

Consider now a reduced-order bilinear model 

* R (k+l) - Aj^x^Ck) ♦ N R [x R (k) • u(k) J + B R u(k) 

y R (k) - C R x R (k) (9) 

where * R (.) I* an n r «l vector, n f < , y R (.) la an n^»l vector, and A^, N R , B R , C R are matrices of appropriate 

dimensions. In addition, we assume that the state covariance X of the reduced model driven by zero ie»a 
Cuasslan white noise Is the unique positive definite solution to 

X R - Vr 4 + (N R°V N R + B R B R • (l0> 

We now define a particular type of reduced order model for discrete bilinear systems. 

Definition : The reduced order model (9), with scate covariance satisfying (10) Is a q- Vol ter ra COVarlance 

Equivalent Realization (q-Volterra C0VER)'of the bilinear system (5) whenever 

\ • V ‘-O - 1 «-> 

and 

R* - R t , 1-0,1 q-1 

where and denote Che sets of order Voi terra and covariance parameters of the 
respectively • * 

An algorithm which constructs the q-Vol terra COVERS of a full order model Is our main 
algorithm is presented next. 


reduced order model. 


objective • One such 
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4 • A Model Reduction Algor ithm 

Suppose Chat a full order nodal (5) and a state covariance satlafylng (6) sre given. The q th obaervablllt; 
matrix ( (3 ] , (4 ) , ( 16] ) of this model Is 


o £ 


^0 


*1-1* 

*1 


*1-1 N 1 

• 

• 

. Qo - c . Qi ■ 

# 

« 


• 

Vi 


*i-l N n 

u 


, i-i <1-1 • 


(u 


The macrlx partitions have dimension (n u +l) 1 'n^*^, 1-0,1 ,...,q-l. We observe that the aatrlcea an 
Q XC*- contain the same Information as the sets and , respectively. Using the full order model w« construe 
t'..e following matrices 


A * 
D - 0 XO 

q q q 


(12 

(13 


o - o (axa**(n[]x)n*)o* - o [ a (NDx l/2 ) )( A (N 0 x 1/2 ) 1 * 0 * . 

q q q q 9 

At a consequence of the quadratic fora and using the Liapunov equation (6) It Immediately follows tbit the rang 
spaces of these matrices are 


R(D ) - R( [ A X 1/2 (NQX I/Z ) B ]) , R(D) - R( [ AX l/2 (NQX 1/Z ) ]) 


(U 


and It Is obvious that R(D ) Is contained In R(D ). 

q q 

We now compute a full rank factorization of 0 


0 - PAP 

q 


(i: 


where rank(D )-r<*n.By virtue of the full rank factorization the columns of P form a basis for the ranj 
q x f 

space of D . Introducing P , the Moore-Penrose Inverse of P, then It la well known that PP la an orthogotu 

projector onto the range of D ([19]). We now partition P Into blocks whose row dimensions are compatible vli 


P - 


q-i 


. p i 


i-l ,...,q-l 


and define new matrices 


q-i 


J-i. 


□ ' ! i 
P q-1 | 


• n u • F " 1 P A P N, 


The matrix G is (n +| )** *n -( n + l)r and tt must be determined such that 
u y u 


P ^ 


D » pXP 


dlag(A) 


n +1 


where I Is a block diagonal matrix with n +1 blocks- Given these constructions we now state our main result. 

u 

Theorem l: Given a discrete bilinear system {A,N,B,C,X} and a matrix G in (17) such that (18) is satisfied tl 

the reduced order model (A.,N ,B C U ,XJ of order n defined by 

R K K K K r 


1 *R N K 1 - * B R 5 P *V • C R ‘ p 0 


X R - A 


c 


where r, P, A are from Che full rank decomposition of D (15), Pg Is from the partition of P (16), and 
satisfies (Id), Is a q-Volterra COVKK. ^ 


Proof: First we will show that P Is the q Ch observability matrix of the reduced order model (19). 


Using 
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decomposition, (IS), (18) and the rang* apaca daaeriptlona (14) wa find Chat 7 la In the range apace of P ao that 
(19) leada to 


P( A R Mr J - P 

which lapllea that Che partition* of P have the required atructure (11) 

P t-1*R 


P 0 • C R • P 1 


P 1-1 N R 


P 1-1 N R 


, 1-1 .... ,q-l 


To ahow that cha reduced order mode] satisfies the bilinear Liapunov equation we first substitute (19) Into 
which leads to 

* e 

A - P*pAP P* 4- P + o BB*0*P + . 

9 q 

Using (12). (13), (13), (18), and by pre and post multiply by P and P , respectively , we have 

* * e 

0 xo • PP 0 AXA 0 P P ♦ PP 0 (NOX)N OP P + PP 0 BB 0 P P . 

qqqq q q qq 

Now ualng the projection property of PP + we find that 

0 (X - AXA* + (NOX)N* + BB*)0* 

q q 

which la known Co be satisfied (6). To show that the model (19) matches Volterra parameters we again use 
projection property 

o B - pp*0 B - OB 

q R R q q 

and the matching of covariance parameters follows directly from (12), (15) 

°„ *a°! ' pf ** - 0 XO* . 

q R R q R q q 

Our remaining Cask Is to determine the unknown matrix C In (17) in order to satisfy (18), This is the 
of the next section. 


5 • Pararaeterlzat Ion of q-Vol terra COVERS 

To obtain a characterization of the matrix G we first examine the structure of the matrices ft and ¥• 
observe that ft ran be partitioned as ^ 

q 

I _ _ 

°q-i d q 

d* 7 

q qq 

and Chat the partitioned form of the conacralnt (18) leads to the three relations 


FlF 


n , , FAG - d , GAG 
q-1 q 


_ * — 
d_ 


qq 


The first relation Is satisfied by virtue of the constuctlon of of P (17). It is easily seen that d 
contained In the range space of F ao that the second relation Is consistent and G* may be expressed as ((19] 

G* - T'CF*! ■*■([- F + F)Y) 

q 

where Y Is an unknown matrix with dimension (n +t)r*(n +1) 1 * *n . Substituting for G In the last relation we 

u u y 

that Y must satisfy the Hermltlan, quadratic, matrix equation 

Y*KY + L*Y ♦ Y*L ♦ M - 0 


X 4 (I - F + F)A 1 ( I - F + F) - K* , L - (I - F + F ) A 


F d , N 

q q 


T 1 f*1 - J 


M 


q qq 

By Inspection we see chat the matrix X Is nonnegative definite, and that the columns of the matrix L 
contal ned in the range space of K. Based on these observations we nnw state a theorem which is motivated 
result of Crone 120]. 


( 10 ) 


Che 

I 

topic 

We 

( 20 ) 

( 21 ) 

Is 

q 

) 

( 22 ) 

find 

( 21 ) 

(24) 

are 
by a 
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Thaoraa 2: Lac K b« an a*a nonnagatlva daflnlta aatrlx with rank t, L an a*n matrix whoaa eoluana ara contalaax 

In Cha ranga apaca of K and N an n*n Haraltlan aatrlx. Than tha aatrlx aquation 

Y*KY ♦ L*Y ♦ Y*L ♦ N • 0 (25: 


haa a aolutton If and only If 

L*K + L - N > 0 and rank(L*K + L - M) • a < t » rank(K) . (24! 

Whan chaaa condition* hold Y la a aolutlon If and only If It has the fora 

Y • K +/2 (VI l/2 U* - K +/2 L) + (I - K +/2 K l/2 )7 (27! 

where K 1 ^ 2 la the unique nonnagatlva daflnlta aquara root of K and K + ^ 2 la tha Moore-Penroee lnvaraa of K*^ 2 . Hu 
aatrlx V la an n«s aatrlx, £ 1a s«a and U la n*s and they auat aatlafy 

V*V - I , R(V) la contained In R(lO , U*U - I , I > 0 , UW* - iVl - M . 

7 la an arbitrary a>n aatrlx. 

Proof : It la wall known chat If K • wfW i j a full rank singular value dacoapoaltlon (SVD), then 

k 1/2 . W(J 1 / 2 U * f k+ /2 . tfn >/ 2 u * 

1/2 ♦ / 2 

and It follows chat K, K , K all have the aaae ranga apace which la spanned by the eoluana of W, an »» 

column unitary aatrlx. By the hypotheals that tha eoluana of L are In the ranga space of K, equation (25) l: 

satisfied If and only If 

(K 1/2 Y + K +/2 L)*(K 1/2 Y + IC* /2 L) - L*K + L - M 
* + 

which la consistent If and only If L K L-M > 0. All aatrlx factors of this relation are 

k i/2 y * k +/2 l - vz l/ 2 u* 

* * 4. 

where Il-.U Is the full rank SVD of L K L-M and V la any colunn unitary aatrlx of appropriate dimension, a «». T 
find a solution Y we must solve Che following linear equation 

K l/2 Y * V£ 1/2 U* - K +/2 L . !2* 

This equation Is consistent If and only If V Is contained In the range spare of K. Since V la coluan unitary :h 
range space of V may he any a dimensional spare with rank a, solutions of (28) exslsc If and only 1 

rank(L K*L-M) * s < t * rank(K). Given that equation (28) Is consistent then Y la a solution If and only If 1 
has Che following fora 

Y ■ K* /2 (VZ ,/2 U* - K +/2 L) ♦ (I - K* /2 K l/2 )Y 


where Y Is an arbitrary ra*n matrix* 


The reaul ts of this theorem show that the matrix L K L-M la the key to solutions of the quadratic macrl 
equation (21). Substituting for X* L, M from equations (24), end using the rules for the Moore-Penrose Inverse o 
a matrix product (121 I), we find that 

lVl - M - 7 - dV T l/2 (I - (F 1/Z (I - F*F))(F l/2 (I - F + F)) + )F l/2 F*d . '29 

qq q q 

From this equation we find an Interesting result on the Moore-Penroso inverse of a quadratic form which we stit 
without proof. 


Fact: The Moore-Penrose Inverse of the quadratic form F~K? is 


(FAfV . F* F l/2 d - (A 1/2 (I - ?*>))() 1/2 d - F + F)) + )A l/ V 


'30 


where T Is a positive definite matrix and F Is any matrix which is niul tlpl Xcatlon compatibles 


Using this result and equation (21) In (29) we find that 


* * — — Jk f — 

L K L - M " d -d D ,d„ M 

qq q q-l q 

which Is guaranteed to be nonnegative definite ((22) )• Thus Che first part of the constraint (2b) of theorem 
wll! always be satisfied. 

To show that the second part of the constraint (2b) will also be satisfied we note that from the range \z*c 
dear rt pt 1 on (14), 
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(32) 


r - rank(Oq) > ranMD^) 

and froa tha definition of K (24), the ralatlona (21) and the roluan diaanalon of F, 

t i rank(K) - rank(I - F*F) • (n ♦Dr - rank(D„ .) . (33) 

Rohde (23] haa shown that tha partitioning (20) of the nonnegative def lnlte matrix lapllaa 

rank(7 ) • rank(D .) ♦ rank(d\ - d*D* ,d ) (34) 

q q-1 qq q q-l q 

and by ualng (31) 

a ^ rankuVl. - M) - rank(0„ ) - rank(0„ .) . (35) 

q q-l 

Collecting aquatlona (32), (33) and (34) we find that a < t and tharefore the aecond part of the conatralnt (26) 
In theorea 2 will alwaya be aatlafled. 

We have ahown that aolutlona of (23), (24) alwaya exlat and by theoraa 2 they will have the fora 

Y - IC +/2 (vr 1/2 U* - K +/2 L) + (I - K* /2 K 1/2 )Y (36) 

where UIU* la the full rank alngular value dcconpoaltlon of L*K*L-M, V la any roluan unitary matrix whoaa range 
apace la contained In the range apace of K and T la arbitrary. We obaerve that tha aecond tern of (36) la In the 

null apace of K which la alao the range apace of F* . It follow* that when (36) la aubatltuted Into the 
expreaalon for G (22) that thla tern will be annihilated by (I-F + F) which repreaenta a projection onto the null 
apace of F along the range space of F* . The first tent of (36) la In the range epace of F, or tna null apace of 
F, so that under the projection (I-F + F) It remains unchanged. Therefore G* becomes 

c* - T'if*! ♦ ic* /2 <vr ,/2 u* - K +/2 U) 

q 

or by ualng equation (24) and conjugate transposing 

* 

c - arV tt'(i - a - f*f)(u - fVjr'u - f^fd^T 1 ) ♦ (37) 

q 

(i - 3V .7 ) 1/2 uv\a - Fair'd - F*F)) +/2 r‘ . 

qq q q-l q 

Equation (37) la an explicit expression for G which was Che objective of thla section. All of the freedom In C 
Is contained In the column unlcary matrix V whose range space la constrained to be In the null apace of F. 

6. Appl lc at Ion to a Robot Mani pu lato r 

Consider the two degree of freedom manipulator Illustrated In Figure 1. The arm has Its center of mass at 
point C, and It may be translated through or rotated about the fixed point 0 by the force F and torque T, 
respectively. The manipulator carries a load at the point L. 



A 

Figure l. Two Degree of Freedoo Manipulator 

Treating the load as a point isass and allowing for Joint stiffness and damping, the equations of aotlon are 
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(■ + M)r" + b f r ♦ k f r - (a + M)r9 2 - Mai 2 • F 

(J ♦ I • Ma 2 ♦ 2Mar * (a ♦ «)r 2 )5' ♦ bgS ♦ k,9 ♦ 2(a ♦ M)rr® ♦ 2Mar« - T 

wh.r. a la tha dlatanca froa C to L, H la cha aaaa of Cha load, J la tha aoaant of lnartla of Cha Joint, a la tha 

aaaa of tha ara and I la lta aoaant of lnartla about C. Joint atlffnaaa and daaping ara rapraaantad by Ic , fc, t a 

h f , bj, raapaetlvaly. Introducing tha atata vactor and tha control 

« • I r r t J | T , u ■ | IT | T 

than tha aquatlona of aotlon have tha generic fora 

* • f(i) ♦ g(i)u . 

A bilinear model of tha manipulator can bo conatructad by expanding each of tha functlona f(.) and g(.) Into a 
power aarlea and Introducing a new atata vactor which containa higher order toraa In t (l*|-|6j). doing the 
flrat threa tana of tha Taylor aarlea axpanalona of f(.) and g(.) and letting 

x • | r r 8 9 r 2 rr 9r 9r r 2 )r Or 0 2 09 8 2 r 2 ... 9 2 ) T 
then wo have a 14 th order bilinear aodel and after dlacretlxatlon it haa tho fora (l). 


For purpoeea of 11 luatratlon, the following mnorlcal valuea are choaen: a * 1 a, a • 100 kg, K • 50 kg, J 
• 1 • 100 kg-a , and k f • 6 N/a, k, - 2.5 N/a, b f • 3 N-aac/a, b, • 3 N-aac/a. Figure 2 ahowa tha atap reaponae 

of the nonlinear equation. of aotlon and tha full order bilinear nodal. Tha bilinear nodal provldea a fair 
approxlaatlon to tha true nonlinear eyatea. A more accurate approximation could bo a.'** by retaining higher 
order terme In the power aeries expansions. 



TIME 

Figure J. Step Response ol Nonlinear and Full Order 3il Inear Models 

Applying the model reduction algorithm with q*3 (matching three sets of Vol terra and covariance parameter*) 
a class of 3-Vol terra COVEKs was obtained. These reduced models have 14 states which Is a greater than fifty 
percent reduction In model order. Figure 3 shows the response of a reduced aodel froa the cJaas of 3-Volterra 
COVERS and the response of the full order model to a unit pulse Input with a 3 second duration. Figure 4 'hows 
the response of thes models driven by a unit Intensity Gaussian white noise process. In Figure 3 we se^ that the 
response of the full and reduced order bilinear models are nearly Identical tor the first 10 seconds. Similarly, 
In Figure 4 the reduced order model mimics the full order model Initially. These observations are In accordance 
with the theory which states chat the response of the reduced order model equals that of the full order systea 
for q steps in time. However, In both cases the quality of the response of the reduced order model deteriorates 
with time and It eventually goes unstable. This Instability Is Input dependent and possibly In a closed loop 
setting the model behavior would be acceptable for greater periods In time. 
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Fi*ur« 1. Deteralnlaclc Reaponae of Full and Kadurad Order Bilinear Model a 
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Figure 4. Stochastic Response of Full end Reduced Order Bilinear Models 


7. Conclusions 


A sequence of sets of Volterra parameters characterizes the detenslnlst if bilinear system, and a sequence of 
sets of covariance parameters describes the stochastic bilinear system. A model reduction technique was developed 
for discrete bilinear systems which generates a class of reduced order models which exactly match the first q 
sets of Volterrs and covariance parameters of the full order model. These models are therefore called q-Volterra 
covariance equivalent realisations, or q-Vol terra COVERS. Methods to choose specific models from within the 
class to satisfy additional modelling considerations is a topic of future research. 
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